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Abstract 



We estimate the di-photon coupling of /o(600), /o(980) and /2(1270) 
i-S^ ' resonances in a coupled channel dispersive approach. The /o(600) di- 

photon coupling is also reinvestigated using a single channel T matrix 
for vrvr scattering with better analyticity property, and it is found to 
(— I ' be significantly smaller than that of a qq state. Especially we also 

estimate the di-photon coupling of the third sheet pole located near 
, KK threshold, denoted as /q^^(980). It is argued that this third sheet 

\ pole may be originated from a coupled channel Breit-Wigner description 

^ ■ of the /o(980) resonance. 

PACS: 11.55.Fv; ll.80.Et; 13.40.-f; 13.75.Lb 

O . 

§ ; 1 Introduction 

The study on the 77 — >■ tttt process has received renewed interests recently on 
^ \ both experimental [1] and theoretical side. [2] One major motivation to drive 

; the studies on the process is to extract the two photon coupling of the resonant 

states appearing in the reaction, which is helpful in exploring the underlying 
structure of these (often weird) states, as emphasized by Pennington. [3] An 
incomplete list of theoretical studies along this line may be found in Refs. [4] 
- [12]. 

In this paper we will also devote to the study on the 77 tt+tt , -k^tx^ 
processes. We proceed with previous calculations, but now take into account 
the coupled channel effects in the most interesting IJ = 00 channel. Section 2 
reviews the coupled channel formalism needed in this study, including: spectral 
function representation, analytic continuation of coupled channel 77 scattering 
amplitudes, and the dispersive integral representation of the latter. In section 3 
we review the basic technique of partial wave expansion and how to estimate 
the Born term contributions. In section 4 we devote to a redetermination to 
the coupled channel tttt, KK scattering T matrices previously proposed, [13] 
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with the help of new experimental inputs at very low energies. [11] Section 5 
is devoted to the numerical analysis on the 77 — > 7r"'"7r^, tt'^tt'^ processes, both 
in coupled channel and in single channel formalism. Finally in section 6 we 
draw conclusions based on the numerical studies. 



2 Couple channel formalism for the 77 — > 
TTTT, KK processes 

2.1 Spectral function representations 

In the coupled channel case the relation between (the 2 by 2) S* matrix and T 
matrix is 

S = l + 2ipi/2ypi/2 ^ 



where p = diag{pi, P2) and p„ = yl — 4m2/s. Here subscripts n = 1,2 
represent vrvr, KK channels, respectively. 

Let Fi, F2 be the 77 — > vrvr, KK amplitudes, respectively, F = (Fi, ^2)"^. 
Define F^{s) = \im^^Q+ F{s ± ie) , then F obeys the following spectral repre- 
sentation in the physical region: 

ImFi = FipiT*, + F2P2T*2 , 

imF2 = f,p,t;, + F2P2t;, (2) 

above KK threshold and 

ImFi = FipiTA , 

ImF2 = F^piT;, (3) 

when s lies between first and second threshold. From Eqs. (2), (3), along 
positive real axis, one has a shorthand expression, 

_ ( K \ _ - ( + 2^Pi^ii^i 22p2Ti2^2 \ ( 



F. 
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- [ 2zpiT2i^i 1 + 2^P2T22^2 [F2 ^ ' 



F^^ = BiiF 



where 9i = 6{s — 4m^) and 62 = 6{s — Amj^) are step functions. 

From above spectral function representations, the analytic continuation of 
F to different sheets can be obtained: 

1 

l+2ipiT/, ' 
l+2ipiTl-^ ' 

^ det 5 ' 
/1 ~2ip2Ti 

F'^ = BjvF = ' i+2^'f^22 ] F . (5) 

V^' l+2ip2T^^ 

These expressions of analytic continuation will be used later when extracting 
the residue couplings of each resonance. 
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2.2 Dispersive representation for 77 — > tttt, KK ampli- 
tudes 

In this subsection we review the method originally proposed by Basdevant 
and collaborators [5] on how to set up the dispersive representation for the 
77 — > TTTT amplitudes. 

On the right hand cut, = SF~. Further let F = F — F^, contains 
only left hand cut of F. In this paper we approximate F^ by its Born term, 
Fb- Then on the positive real axis one has, 

F = F + Fb , (6) 
F+ = {F + FbY = F+ + Fb = SF- = S{F- + Fb) , (7) 
F+ = SF~ + {S- 1)Fb . (8) 

Next we will search for a 2 x 2 invertible matrix function D{s), which is analytic 
on the cut s plane and only contains right hand cut {r.h.c), and on the r.h.c. 
it obeys the same unitarity equation as F{s). That is 

D+{s) = S{s)D-{s) , or S = D+{D-)-^ . (9) 

One can then deduce from the above equation and Eq. (8) that 

ImD-^F = -IirD-^Fb , (10) 

from which one obtains a dispersive representation for F~^ 

F+ = -— "" ds' + D+P, (11) 

TT Jso s' — s — ie 

where P is a 2 dimensional array of subtraction polynomial, and a proper 
subtraction is understood on the dispersion integral. By making use of Low's 
theorem [16] which tells that when s — 0, F{s) — > Fb{s), we can rewrite the 
above equation with one more subtraction, and obtain: [17] 

F(s) = + D(s)[Ps - f! r ■°;f;<-'>^^<-'' d/] , (12) 

From the above discussion we know that the 77 ^ vrvr, KK amplitudes can be 
fit with parameter (s) P, once D{s) and Fb{s) are known. 



2.3 Solution of function D{s) 

The 2 by 2 matrix function D{s) only contains right hand cuts, and satisfies 
the similar unitarity relations as F: 

ImDn = D,^p,Tl^d^ + D21P2TI2O2 , 

lmD2i = DiipiT;,9i + D2iP2T;^e2 ; (13) 
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lmD2i = DupiT^iOi + D22P2TI2O2 , 

ImD22 = Di2piT;,ei + D22P2T22O2 . (14) 

The above two equations have the same structure as the following 

Imxi = XipiTliOi + X2P2TI2O2 , 

ImX2 = X1P1T21G1 + X2P2T22O2 ■ (15) 

Hence searching for solutions oi D{s) is equivalent to searching for two indepen- 
dent solutions of the equation for the two dimensional array (xi)X2)^- The 
two independent solutions can be identified as {Du, D2i)^ and (Di2,-D22)^- 
The integral equation to solve Eq. (15) is 



IT J Ami S'{s' — S — le) 

s f^^ Re[x^{s')p2{s')T*,{s')92]^^> 



TT J 4ml S [S — S — le) 

TT J Ami s'{s' — s — ie) 
^ s rA^ Re[x^{s')p2{s')T;,{s')e2] ^, ^ 
TT Jiml s'{s' - s - ie) 

where denotes the number of steps in the iteration, the integration is trun- 
cated at = 2.4GeV^. On the r.h.s. of above equation, we take the real 
part of the numerator in the integrand which is of great help in increasing the 
speed of convergence. [18] In the numerical calculation a convergent solution 
emerges after approximately 15 steps of iteration. 

In this paper we solve the coupled channel D{s) function in the case of 
/, J = 0, 0, for a given T matrix. For the two d waves and the I = 2 s wave we 
only use single channel approximation, and in such simplified case the solution 
to the function D{s) is known as the Omnes solution, 

2.4 Analytic continuation, pole residues and r(/ ^ 77) 

In order to extract pole residues on different sheets a knowledge of analytic 
continuation to the complex plane is needed. Firstly, the analytic continuation 
of function D{s) is simple. It satisfies the following dispersive integral equation 

TT JAml S'[S' - S) 

D{s), p{s),T{s),6 in above can all be matrix functions in above formula. To 
calculate D with complex argument z on the first Riemann sheet one only 
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needs to replace s ^ z in above integral formula. For function F the analytic 
continuation is more involved, in the following we discuss this topic in some 
details. 

Residues of second sheet poles: 

Since in the vicinity of a pole zji, S{s ~ Zfj) ^ S'{zii){s — zjj) and 

T^^inn n7r){s - Zjj) = S^\7'''\ (19) 

S'{zii){s - zn) 

F^^(77 ^^n){s^zn)= J/^\T '''\ (20) 

S'{zii){s - Zn) 

and through the definition of coupling constants: 

^2 

T(7r7r ^ 7r7r)(s — s> 2fl) ^ 



Zr- s ' 



F(^^-^Tr7r){s^zn) = , (21) 

zr- s 

one obtains the expressions of and g^y for second sheet poles: 
2 Tjzu) Fjzii) 
S'{zn) S'{zn) 
Residues of third sheet poles: 

Denote the third sheet pole as zm, which is the zero of det 5, then 

det S{s ^ ziii) = (det S)\ziii){z - zui) ■ (23) 

From 

rflll _ T^III _ M 11 ^12 det a, Acts 1 (')A\ 

the residues are hence obtainable, according to Eq. (21): 

2 ^ 'S'22(^m) 2 ^ 'S'll (^JJj) ,c2^rC\ 

2zpidet^'(27//) ' 2ip2det5'(^///) ' ^ ^ 

From Eq. (5) one further gets 

Fl{zin)S22{ziii) Fl{zHi)2ip2{zni)Tl^{zni) , . 

^"^^ (det 5)'(^,,,) + (det 5)'(z,,,) ' ^'^^ 

^ Fj {ziii)2ipi {zni)Tl2 (zju) _ F^{zjjj)Sn{zjii) 

^''^"^ (det Syizrn) {det Sy{zni) ' ^ ^ 

Having obtained the two photon couplings, g-y, the width of a, /o(980) and 
/2(1270) to two photons can be calculated, 

r« (pole) = "'^^'^^1' (28) 
^^y[Po^e) 4^2J + l)m^- 

where a = 1/137 is the fine structure constant and /3r = {1 — 4m^/m^)^/^. 
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3 Partial Waves and the Born term contribu- 
tions 



3.1 Partial wave expansions 

For the 77 — > tttt process there are two independent hehcity amphtudes M_|.+ , 
M+_. The former corresponds to photon hehcity difference A = the latter 
corresponds to A = 2. They contribute to the differential cross-section as: [ 1] 

where p = [1 — 4m^/s]^/^, and they have partial wave expansions involving 
only even J(> A), 



M+4s,e,(j)) = e'Vl67rJ2 Fjo{s)Yjo{e,(j)) , (30) 

,/>o 



M+_(s,e,</.) = e2v^^Fj2(s)y;2(^,0) . (31) 

J>2 

With this normalization the total cross-section reads, 

^ ,j>x 

At very low energies the major contribution to the amplitudes comes from 
one 71 exchange (OPE), i.e., the Born term, as guaranteed by Low's theo- 
rem. [16] 

M+±(s,t)^M^i(s,t), {s^O,t^ml). (33) 

Denoting s-channel partial wave of the Born amplitude as Bjx{s), then Low's 
theorem implies 

^^-1, Fj,{s) - Bj,{s) 0{s) . (34) 
Bjx{s} 

The OPE Born term amplitude for 77 tt^tt^ vanishes and for 77 ti^ti~ 
processes read, 

= = 2eW(^^ + , (35) 



„ R 26^ {m^ — tu) , 1 1 , , ^ 

= Mf^ = ^ + — ) ' (36) 

s — t m'^ — u 

from which OPE partial wave amplitudes are obtained: [1] 

BooW = ^^ln(ii^), (37) 

j(l^|P^,„(i±^)-6l. (38) 

V lo p^ p 1 ~ P 

^ , , /T5il -p2)2 , ,1 +p, 10 2 ^ 

^22 ^ = \ —[- r^ln ^ -) + ^ - - • 39 

V 32 p'^ 1 — P 3 
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At higher energies, however, OPE Born term amphtudes are not satisfac- 
tory to describe the left cut well. One may try to improve this by adding 
crossed channel vector and axial vector meson exchange diagrams. Never- 
theless it is difficult to judge to what extent adding vector and axial vector 
exchange contributions improve the situation. Further discussions will be made 
in section 5.1. 

4 The coupled channel scattering T matrices 

It is necessary to choose an appropriate T matrix in each channel. For the 
1=2 s-wave we use the result of Ref. [19]. For the 1=0 d-wave amplitude we 
use the result from Ref. [20]. In both these two cases we use single channel 
approximation. 

For the 1=0 s-wave vrvr ,KK system, to solve the matrix function D{s) 
by iteration, a general coupled channel T matrix is required as an input. We 
choose the i^- matrix parametrization proposed in Ref. [13] (the form) to 
refit low-energy vrvr scattering data with totally 22 free parameters. The data 
sets we use are the phase shift and inelasticity of / = s-wave vrvr scattering 
from CERN-Munich [21] and especially from the K^i experiment [11],^ and 
012 = + of HIT — >• KK 1 = s-wave scattering from several groups. [22, 
23, 24, 25, 26] We have total 171 data points below 1.89GeV in the region of 
our interest. It is worthy mentioning that we can get an acceptable fit with 
a Xdof 1.14, even including the confiicting data of 0i2 from Cohen et al. 
and Etkin et al., but the exercise leaving aside either of them favors Etkin et 
al. over Cohen et al. with a. Xdof 0-68 compared with 0.99. The final fit 
excluding the data from Cohen et al. are shown in figure 1 and figure 2. The 
obtained i^-matrix parameters are listed in table 1. 




Figure 1: The fit curve of vrvr J = S-wave phase shift and inelasticity with 
CERN-Munich data [21] and data from Pislak et al. [14]. 



^Adding the data from Ref. [14] is already good enough to constrain the T matrix near 
threshold so the most recent NA48/2 Ke4 data is not included. [\J>] 
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Figure 2: The phase shift 0i2 = + 5k of vrvr — ^ KK 1 = s-wave scattering 
with data sets from Res. [22, 23, 24, 25, 26]. Notice that Ref. [22] is not used 
in the fit. 



So 


= 0.081 


4i 


= 7.599 


'-21 


= 15.49 


Si 


= 0.305 




= -16.30 


'^21 


= -15.60 


S2 


= 0.984 


'-11 


= 7.530 


'-21 


= 17.58 


S3 


= 3.587 


'-11 


= -2.796 


„3 
'-21 


= -4.942 


n 


= 1.499 


'-11 


= 0.0 


'-21 


= 0.0 


f? 


= 1.041 


„o 

'-12 


= 15.49 


„0 
'-22 


= -6.441 


n 


= 0.089 


'^12 


= -15.60 


'-22 


= -13.53 


/I 


= 0.484 


'-12 


= 17.58 


'-22 


= 31.04 


/I 


= -7.714 


„3 
'-12 


= -4.942 


„3 
'-22 


= -6.830 


/I 


= 9.216 


'-12 


= 0.0 


^-22 


= 0.0 



Table 1: matrix parameters obtained from the fit as described in this 
section. For the definition of these parameters we refer to the original paper 
Ref. [13]. 



The fit to the amplitude produces several S-matrix poles either on the 
second Riemann sheet or third-sheet. Those that are relevant to the current 
interest are listed in table 2. There we also list a third sheet pole located at 
0.705 — 0.3272, which might be considered as the third sheet counterpart of 
/o(600). Since it is too far away from physical region, we will not discuss it any 
further. The fit with a i^s parametrization returns another pole on the third 
sheet, denoted as /q^^(980) herewith. It is of course a {K^) model dependent 
prediction. This third sheet pole is close to but below KK threshold, and may 
be correlated to /o"(980), [13, 27, 28] since the data are better described when 
including such a pole.^ We will further discuss physics related to this pole in 
some detail in section 5.3. 

^It is worthy emphasizing that the third-sheet pole near KK threshold moves away and 
disappears, if the data of Cohen et al. are taken into account. 
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pole 


sheet-II 


sheet-Ill 


a 


0.549 - 0.230Z 


0.705 - 0.327i 


/o(980) 


0.999 - 0.021Z 


0.977 -O.OeOz i 



Table 2: The poles's location on the ^/s-p\a.ne, in units of GeV. 



pole position 


JTTTV 


9kk 


y/sj] = 0.999 - 0.021i 


-0.07 -O.OH 


-0.10 + 0.092 


^sin = 0.977 - O.mOt 


-0.10 + 0.02i 


-0.02 - 0.092 



Table 3: /o(980) pole positions and their residues. 



The a pole is not very satisfactory comparing with the determination of 
Ref. [19], since a coupled channel vrvr, KK scattering amplitude fully compati- 
ble with analyticity and crossing symmetry is still un-available. In section 5.2 
we will try to remedy this shortcoming by making use of the single channel T 
matrix of Ref. [19]. The pole position is not very satisfactory, and the coupling 
strength extracted, 

gl^ = i-Om-0.m)GeV\ (40) 

differs from that given in Eq. (43) in section 5.2, though the magnitude is 
compatible. On the other side it is expected that the information one extracts 
from the coupled channel is reliable in the vicinity of KK threshold, especially 
for the /o(980) resonance. Couplings to vrvr and KK of poles near KK thresh- 
old are also obtained as listed in table 3. It will be useful in section 5.3 when 
discussing the properties of /o(980). 

5 Numerical fit to 77 — > tt^tt", tt'^tt'^ process 
5.1 The coupled channel fit 

Only I,J=0,0 channel are calculated by solving coupled channel integral equa- 
tions. Other channels are all approximated by Omnes solutions. Subtraction 
polynomial (constant) in different channels are listed below: 

• I, J, X = 0, 0, 0: Solve couple-channel integral equation for D and fit 
parameters are Pi and P2. 

• I, J, \ = 2, 0, 0: Use Omnes solution, the fit parameter is P3. 

• /, J, A = 0, 2, 2: Use Omnes function and fit P4. 

• /, J, A = 2, 2, 2: Since this channel's contribution is very weak, we use 
Born term approximation. 

• /, J, A = 0, 2, 0: Use Omnes solution and the fit parameter is P5. 
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• /, J, A = 2, 2, 0: Since this channel's contribution is very weak, we use 
Born term approximation. 

Here five subtraction constants are involved. Furthermore we introduce two 
additional form-factors to suppress the bad high energy behavior of the Born 
term amplitudes: 

Bs^B.e-''^', 5,^5,e-/A^, (41) 

hence totally we have 7 fit parameters. 

For 77 — s> TT+TT" process there exist three sets of data, which are from 
Belle, [1] CELLO, [29] and Mark-II. [30] In our fit from Tr+Tr" threshold to 
0.9GeV, we use the Mark-II data, from O.DGeV to 1.4GeV we use Belle data. 
For 77 TX^-K^ there exist two data sets from Crystal Ball Collaboration, 
Refs. [31] and [32]. Here we chose the data from Ref. [31], from tttt threshold 
up to 1.4GeV. 

Firstly a comparison between the contribution from OPE Born term ampli- 
tude and the contribution adding vector and axial vector meson exchanges in 
each channel is made, having obtained the two cutoff parameters which used 
to readjust the Born term amplitudes: = 1.77GeV and = 0.86GeV. 
However it is difficult to judge wether adding vector meson and axial vector 
meson exchange contributions improves the OPE contribution or not, com- 
paring with the fit results shown in figure 3, where contributions from one tt 
exchange, vector meson exchange, axial-vector meson exchange, with correc- 
tions from the exponential form-factors (the latter are determined from the 
fit), are combined together and plotted. For the expressions of vector and 
axial vector meson exchange contributions we refer to the appendix. 

The ratio between IJ\ = 020 wave and IJX = 022 wave may not be very 
stable. Though it is understood that the former should be much smaller than 
the latter. [33] We however find a convergent solution with (Jdo/<^tot — 0.14 and 
(^slc^tot — 0.07 (at the /2(1270) pole position). The first ratio given here is 
rather close to the 'dip'' solution of Ref. [31]. The second ratio is significantly 
smaller than both the 'peak' and the 'dip' solution of Ref. [31]. The fit results 
are listed in table 4 and plotted in figure 4. The fit result on each partial 
wave amplitude is found to be rather similar in general with the solution B of 
Ref. [2] and to figure 1 of Ref. [9]. The di-photon width r(/2 27) = 4.36keV 
here, which is larger than the value given by the solution B of Ref. [2]: r(/2 — 
27) = 3.82 ± 0.30keV. Part of the reason for such a difference may be due to 
the fact that the s-wave contribution given in this paper at the /2(1270) peak 
is smaller than that given in Ref. [2]. Also the interference between s wave 
and d wave is found to be destructive.^ That the result of Ref. [2] on /2(1270) 
di-photon width is significantly larger than the value as quoted by PDG is 
explained in Ref. [2] - because the Belle data gives a larger enhancement at 
/2(1270) peak than the previous data. 

■'Due to the limitation of experimental detection, the integration range of cosO is not 
from -1 to +1, hence causing a nonvanishing interference between different partial waves. 
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Figure 3: Born term amplitudes with corrections from exponential form factors 
taken into account. The parameter of form- factors are determined from the fit 
made in section 5.1. 



Pole-positions(GeV) r(/j ^ 77)(keV) 



^P(980) 0.999 - 0.021Z 012 

/o^^^(980) 0.977 - O.OeOz 0.35 

/o(600) 0.549 - 0.230Z 0.76 

/2(1270)(A = 0) 1.272 - 0.087i 0.66 

/2(1270)(A = 2) 3.70 



Table 4: Xdof~ 0-^' using T matrices as depicted in section 2.3. The central 
value of other 5 parameters are: Pi = 0.25, Pa = 0.25, P3 = -0.23, P4 = 0.22, 
P5 = 0.55. 



To summarize, we have the following observations through our fit: 

1. It is known that the A = 2 amplitude should dominate the 77 width of 
the tensor state. Even though the ratio between the A = and the A = 2 
partial wave may not be stable, a solution with a small ratio ado/o'tot is 
found. The total di-photon decay width of /a (1270) is however found to 
be rather large comparing with the value quoted by PDG and the value 
given in Ref. [2]. 

2. The /o(980) — > 77 width is found to be very small. It is in agreement 
with the solution B in Ref. [2]. See table 6 for a comparison between 
different model predictions. 

3. We also give the /o^^(980) width to two photons, which is not seen in 
previous literature. However the third sheet pole's position and hence 
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Figure 4: The coupled channel fit to the 77 — * tt+tt , tt^tt^ data. The data 
sets chosen are depicted in the text. 

its residue are not quite stable, hence this number should be treated with 
care. 

4. The (J width is smaller than most values given in the literature. 

5.2 A refined analysis on the sigma coupling 

We have noticed that the low energy physics related to the a pole as described 
by the coupled channel T matrix is not very satisfactory. Because the lack 
of crossing symmetry, occurrence of spurious poles, and the distortion of the 
o pole location. All these defects could affect a reliable extraction of the two 
photon coupling of the a meson. We try to remedy these defects by refitting the 
low energy data using a single channel T matrix with all the nice properties 
such as analyticity, crossing symmetry, unitarity and with a reliable a pole 
location, since there is no coupled channel T matrix with these properties 
available yet. The appropriate choice is the tttt scattering T matrix proposed 
in Ref. [19], which gives the a pole location (The Eq. (21) of Ref. [19]): 

= 457MeV , = 551MeV , (42) 

and the residue: 

gl^ = i-0.20-0.m)GeV\ (43) 
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Xci.o. f. 


Pi 


Pole-position (GeV) 


T{a- 


77) (keV) 


0.8 


0.49 


0.457 - 0.276Z 


2.08 



Table 5: A single channel fit up to 0.8GeV data using T matrix of Ref. [19], 
with only one parameter Pi. 



to be compared with Eq. (40).^ Also the value in Eq. (43) is compatible with 
the result of Ref. [12, 36], gl^^ = -0.25 - OMiGeV^. 

We refit the 77 Tr+Tr" , 7r°7r° data below 800MeV using that Tj=^ and 
Tj^Q matrices and keeping d-w&ves fixed by the coupled channel fit described 
in section 5. That means contributions from all other partial waves are treated 
as background. In this situation we fit with only single parameter Pi, while 
keeping all other parameters fixed by the coupled channel fit. The result is 
given in table 5. We plot the fit curve in figure 5. 




Figure 5: A fit up to 0.8GeV using single channel s-wave T matrices of Ref. [19] , 
with one fit parameter Pi. The tt+tt^ and tt^tt^ data are from Refs. [30] and 
[31], respectively. Dashed curve represents d-wave background, solid curve 
represents the total contributions, including the 1=0 s-wave to be fitted. 



5.3 A possible Breit— Wigner description to the /o(980) 
resonance 

In this section we discuss the possibility that the two narrow width poles in the 
vicinity of KK threshold found in section 4 may be originated from a single 
Breit- Wigner parametrization. 

^Thc g^TTTT coupling differs by a factor 16n from that of Ref. [Id] where the gaTnr — 3GcV. 
Ref. [1(1] also quoted the number extracted from Ref. [ ■ ] in two ways, corresponding to 
IgaTTirl = 3.4 and 3.9GeV, respectively. 
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For the vrvr, KK coupled channel system, one can make use of conformal 
mapping technique [37] to map the four sheets s plane into the one sheet u 
plane, 



(44) 



where pi = s — 4m? denotes the i-th. channel momentum, A = ^7722 — mf. 
The cu-plane is depicted in figure 6: In figure 6 uj = i,l corresponds to s = 4m^, 





i 


I 




II \ 


1 


1 


IV 1 


III 



Figure 6: Conformal mapping and cj-plane. 
4m|f, respectively. S matrix elements can be written as. 



d{uj) 



det S{uj) 



d{-uj) 



(45) 



Function d{u) contains no kinematical cut and d{—u!*) = d*{uj), according to 
real analyticity. An S matrix pole corresponds to a zero of d{uj). We denote 
second sheet pole and third sheet pole by uj2 and UJ3, respectively and 



UJ2 
Us 



r2exp [i(p2) 
rg exp (i03) 



ra < 1 , 02 > 
rg > 1 , 03 < . 



The relations between r^, 0j and s-plane parameters S2 = (M2 
(Ms - |r3)2 are 

MiTi -- 
M^-Tl/A - 
In such a two-pole case. 



(46) 

S3 = 



(-^-r/)A^sin: 



(— + r2)A2 cos: 



+ 2{ml + ml) 



where 



d{u) = d^3[u)D{uj) 



(47) 

(48) 
(49) 
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and then 

= — dR;^ = ^(Z^T ^ ^ ■ ^ ^ 

One can prove, under some simplifications, that when r2r3 = 1 the two 
pole description Eq. (50) is equivalent to the following coupled channel Breit- 
Wigner description [38], 



So - s - i{-pi^ + p2^ 



So - s- i{pi^ + P2y) 



exp {2i5bg) , (51) 



where 



So = 4m^ + A2[(r3- — )2 + 4sin02sin03] , 
7,r = 4m2A(r3 - — )(sin02 - sin^s) , 

7fc = 4m2A(r3 + — )(-sin</)2-sin03) . (52) 

The second equality in Eq. (51) hold approximately, because we limit our- 
selves in the energy region around KK threshold. Using the pole positions 
provided by table 2 and Eq. (47) we find the pole locations on the oj plane as 
following: 

ra = 0.82, sin02 = O.12, 

rg = 1.27, sin03 = -O.28, (53) 

from which one finds r2r3 = 1.04 ^ 1. This suggests that the table 2 does imply 
that /d^(980) and /o^^(980) may come from the same Breit-Wigner resonance. 
The Breit-Wigner resonance parameters in Eq. (51) can be determined (taking 
l/r2 as r^) using Eq. (52), 

So = 0.982 GeV^ 7^ = 0.146 GeV^ 7j^ = 0.304 GeV^ . (54) 

Moreover, using the information obtained from the T matrix fit on /q ^ coupling 
to TiTX (see section 4), 

[gl)jj = (-0.071 - O.Oll'OGeV^ , (55) 

and comparing it with Eq. (51) one can actually determine the background 
phase Sbg{s) at = 2mK to be ~ 91.8°,^ which is actually mainly contributed 
by the a pole. [19] In this determination one make use of the fact that fo is 



^This background phase has also been discussed in Ref. 
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a narrow width resonance and hence one can assume Sbg{s) calculated at the 
pole position is approximately real. 

Table 4 further indicates that the /q^^'s di-photon coupling is much larger 
than that of /q^: |(5'^)m/(5'^)//| — 3. On the other side an analysis based on 
coupled channel Breit-Wigner description (Eq. (51)) and the two pole positions 
provided in table 4 gives the ratio to be ~ 1.5. Though the two ratios are rather 
different, nevertheless they at least both predict \ {g^)ni\ > \{g^)ii\- Hence the 
two narrow width poles in the vicinity of KK threshold found in section 4 may 
still be considered as originated from a single Breit-Wigner parametrization, 
in support of an early suggestion of Ref . [28] . Of course much more efforts still 
need to be done to clarify such a issue. 

6 Discussions and conclusions 

A coupled channel dispersive analysis on the 77 tt+tt", Tr^TT*^ processes is 
made using a K matrix parametrization from Ref. [13], but refitted with the 
new Ke4 data of Ref. [11]. The shape of different 77 — vrvr partial wave cross- 
sections are similar to those given in Refs. [2, 9], though we have a smaller s 
wave contribution at /2(1270) peak comparing with Ref. [2]. Properties of two 
poles (one on sheet II, one on sheet III) found near KK threshold are investi- 
gated and it is found that the two poles may be explained as coming from a 
single Breit-Wigner parametrization. Our prediction on the two photon width 
of the second sheet /o(980) resonance as listed in table 4 is smaller comparing 
with most previous determinations found in the literature, but agrees with the 
solution B of Ref. [2]. 

A refined analysis on the di-photon coupling of /o(600) is made, using the 
single channel T-matrix (the PKU) parametrization established in Ref. [19, 
40]. The PKU parametrization maintains the nice property such as unitarity, 
analyticity and (numerically) crossing symmetry. Our result in table 5 gives 
r((T — > 27) ~ 2.1keV, which is compatible with the result of Ref. [10], though 
in latter the investigation only is only confined in the energy region below 
O.SGeV and the d-wave contribution is not considered. 

However, by comparing table 4 and 5, we find that the di-photon coupling 
of the a pole may not be very stable. Nevertheless one may argue, based on the 
analysis made in this paper, that the coupling ought to be significantly smaller 
than the estimate based on a simple {uu + dd)/^/2 assignment. We borrow 
from Ref. [48] table 6 on radiative width of scalars in different modelings of 
their composition, adding the result from Ref. [43]. Clearly one reads from 
table 6 that a scalar with di-photon decay width significantly smaller than 
4keV cannot be of a simple qq nature.^ It is argued in Refs. [50, 51] that the 
/o(600) is the chiral partner of the pseudo-goldstone boson, or the a meson 
responsible for spontaneous chiral symmetry breaking. In this picture one 
expects the /o(600) meson can not be described simply by a pure qq or a 

^See however Ref. [49]. 



16 



composition 


prediction 


author (s) 


(uu + dd)/y/2 


4.0 


Babcock and Rosner [11] 


ss 


0.2 


Barnes [12] 


99 


0.2 ~ 0.6 


Narison [13] 


[ns] [ns] 


0.27 


Achasov et al [11] 




0.6 


Barnes [45] 


KK 


0.22 


Hanhart et al [46] 



Table 6: Summary of two photon decay width of scalars calculated in different 
models. 



pure q^q^ picture, rather it is a mixture of many components of Fock state 
expansion, possibly includes a sizable glue content as well. [12, 36] It is nature 
then to expect that this property is reflected by its two photon coupling as 
suggested by the estimation made in this paper. 
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8 Appendix 

In this section we describe the method how we calculate various vector meson 
and axial vector meson exchange contributions to the Born term amplitude. 
We use Proca field to describe spin 1 fields. For the 1 vector meson, the 
interaction lagrangian is 



kin 



1 

4 



(56) 



Proca field is also suitable to describe the two types of axial vector meson 
contributions: A^iJ^^ = 1++) and ^^(J^^ = !+")• 

For TTirj vertices, the lowest order amplitudes are calculable using the 
following interaction lagrangian: [52] 



C^nt{V) = eCve^^p^F^" <V,{Q,u''}>, 



(57) 



where 



V,.. 



A,. 



B,. 



/ Po I ^ 



V 



p 

K*- 



/ , /i(1285) 

V2^ V2 



k: 



lA 

/ , fei(1170) 

hi 



'IS 



K*+ \ 

-0 ; 



V2 + 



/i(1285) 



bt 



V2 



b'l fei(1170) 



^lA 
^lA 



mm), 
Kts ^ 



/ii(1380) j 



(58) 



The KiA and Kib mesons are mixed states. 



KiA = sin e/sTi (1270) + cos ^A'l (1400), 
KiB = -sin^iri(1400) + cos^Ki(1270). 



(59) 
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The mixing angle |^| = 37° or 58° [53, 54, 55]. Comparing with Ref. [52], we 
also introduced exchanges here except the 1 , l'^" contributions. 

Coefficients Cy ,Ca,Cb are determined through p(770) 717, ai(1260) n'j 
and 61(1235) — > 7r7 processes. 

Invariant amplitudes A, B are obtainable using the above interaction la- 
grangian: 

•77 ~^ 7r°7r° : pQ,uj,h\,hi{inQ) contributes here. 

9F2 Mj-t ' 



B 



t 



^7v,<^ - — -t — ~ ^^nA^^-^p ^ ^^-^^f' 
^NM 35^2 Ml - 1 ' 



72F2M2 -t' 



^k/^i = ^]^^^ = 9^kfti (^^1^^/^1(1170)), (60) 



To get u channel contributions we only need to replace t by m in above expres- 
sions, i.e. Aj^^i = A%^i{t u), B^i = B%^^{t u) 

•77 — i> TT+TT^ : vr, p, 61, ai contribute here. 

A Contact ^ 



TjContact n 



A: 



{2t + s - 2m: 



2\2 



t 



^2 ^^2 _ ' 



5^ 



^ C'y s-At- Ami ,t 

^C,p ^g^2 ^2 _ ^ ^A^,P' 
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R< — B — R* 

^c,ai = -^^^^^ = -9i?^,,,(CB^CA,M,, ^M,J, (61) 

To get u channel contribution we only need to replace t by u in above expres- 
sions, i.e. = A'cA^ n), B^^^ = B^^,{t ^ u) 

•77 K+R- : K,K*, Kib, Kia contribute. 
4 



A Contact 
^C,K 

TjContact r\ 
^C,K = U, 



s 



, _ (2t + g-2m^)^ 1 _ 
^c,i^ = —2—: = Bc^^m^ uik), 

1 

(^2 ^ _|_ AuP 

^c,K,, = 3g^2 ^2 _ = ^ m^,M,, ^ Mk,„F^ Fk), 



^c,/fi, = i^J^T^—^=B%^,,{m^^mK,M,,^MK,,,F-.FK) 



(^2 ^ _|_ 4^ _ 4j^2 



(62) 



To get M channel contributions we only need to replace t by u in above expres- 
sions, i.e. A^c, = ^cAt ^ u), B^. = BUt ^ u) 



77 K^K : JC, KiB contribute. 

4(^2 _g _ 4^ _ 4^2 

= MK-t"" ^ ^^nM- ^ ^i^. M, ^ M^., F ^ Fk) 

= Ml^-t ^ 4E^.p(^^ ^ ^i^^ Mp^Mk*,F-^ Fk), 
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(63) 

To get M channel contributions we only need to replace t by n in above expres- 
sions, i.e. = Al,.{t u), B^^ = B'c^t ^ «)• 

In above formulas, in a term ^^"^^j, the subscript N{C) denote neutral 
(charged) particle exchanges. Superscript i means resonance Ri contribution 
in such channel, t{u) denotes t(M)channel. 

One can hence get helicity amplitudes H^^ for processes 77 — > 7r"'"7r~, 
77 ^ vtOttO, 77 ^ /ST+K", 77 ^ fsfOZ": 

= A + 2(4m2-s)5, 

H^. = ^J<^B. (64) 
s 

For example, for the 77 tt+tt^ process, the helicity amplitude read, 

TT /) Contact I A t+lL 1 A t+U I ^ t+« I A t+u 



— tu) 



^ -V--. (ijg-^-^ + + + B'+l + ) . (65) 

The amplitudes are expressed using Condon- Shortly convention. 

M+_ used in this paper are defined using non Condon-Shortly convention 

2 2 
M++ = , M+_ = ^^H^- . (66) 

Finally we plot in figure 7 various Born term contributions in ti^ti~ and 
TT^TT^ channels. 
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